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Abstract 

Background: The mitogen-activated protein kinase (MARK) cascade is an evolutionarily ancient mechanism of 
signal transduction found in eukaryotic cells. In plants, MARK cascades are associated with responses to various 
abiotic and biotic stresses such as plant pathogens. MARK cascades function through sequential phosphorylation: 
MARK kinase kinases (MARKKKs) phosphorylate MARK kinases (MARKKs), and phosphorylated MARKKs phosphorylate 
MARKs. Of these three types of kinase, the MARKKKs exhibit the most divergence in the plant genome. Their great 
diversity is assumed to allow MARKKKs to regulate many specific signaling pathways in plants despite the relatively 
limited number of MARKKs and MARKs. Although some plant MARKKKs, including the MARKKKa of Nicotiana 
bentliamiana (NbMARKKKa), are known to play crucial roles in plant defense responses, the functional relationship 
among MARKKK genes is poorly understood. Here, we performed a comparative functional analysis of MARKKKs to 
investigate the signaling pathway leading to the defense response. 

Results: We cloned three novel MARKKK genes from N. bentliamiana: NbMAPKKKP, NbMAPKKKy, and NbMAPKKKe2. 
Transient overexpression of full-length NbMARKKK(3 or NbMARKKKy or their kinase domains in N. bentliamiana 
leaves induced hypersensitive response (HR)-like cell death associated with hydrogen peroxide production. This 
activity was dependent on the kinase activity of the overexpressed MARKKK. In addition, virus-induced silencing of 
NbMAPKKK(i or NbMAPKKKy expression significantly suppressed the induction of programmed cell death (RCD) by 
viral infection. Furthermore, in epistasis analysis of the functional relationships among NbMARKKK(3, NbMARKKKy, 
and NbMARKKKa (previously shown to be involved in plant defense responses) conducted by combining transient 
overexpression analysis and virus-induced gene silencing, silencing of NbMAPKKKa suppressed cell death induced 
by the overexpression of the NbMARKKK(3 kinase domain or of NbMARKKKy, but silencing of NbMAPKKK(5 failed to 
suppress cell death induced by the overexpression of NbMARKKKa or NbMARKKKy. Silencing of NbMAPKKKy 
suppressed cell death induced by the NbMARKKK(3 kinase domain but not that induced by NbMAPKKKa. 

Conclusions: These results demonstrate that in addition to NbMARKKKa, NbMARKKKP and NbMARKKKy also 
function as positive regulators of RCD. Furthermore, these three MARKKKs form a linear signaling pathway leading 
to RCD; this pathway proceeds from NbMARKKK(3 to NbMARKKKy to NbMARKKKa. 
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Background 

Because plants lack an adaptive immune system, appro- 
priate perceptions and responses of individual cells to 
various environmental stimuli, such as the biotic stress 
caused by phytopathogenic microorganisms, are critic- 
ally important. The plant defense response against biotic 
stress is triggered by the recognition of conserved 
pathogen-associated molecular patterns (PAMPs) or of 
pathogen strain-specific factors known as elicitors or 
effectors [1]. The response triggered by PAMPs is known 
as the basal defense response, whereas that triggered by 
specific elicitors is known as the hypersensitive response 
(HR). In the latter, an effector is recognized by a corre- 
sponding plant resistance (R) protein. The HR is fre- 
quently accompanied by programmed cell death (PCD), 
which plays a particularly important role in the defense 
against biotrophic pathogens but is also an essential 
function in normal plant development and differenti- 
ation [2]. Although many plant components required for 
the PCD-associated HR have been identified, the entire 
signaling pathway leading to PCD has not been elucidated. 

The mitogen-activated protein kinase (MAPK) cascade 
is a highly evolutionarily conserved signal transduction 
mechanism found in eukaryotic cells. Subsequent to ac- 
tivation of the cascade by various extracellular stimuli, 
the signal is transduced intracellularly by sequential 
phosphorylation. In plants, MAPK cascades are asso- 
ciated with developmental and hormonal responses and 
with stress responses to abiotic and biotic factors [3]. A 
MAPK cascade consists of three functionally linked pro- 
tein kinases: a MAPK is phosphorylated and activated by 
a MAPK kinase (MAPKK), which is in turn activated by 
an upstream MAPK kinase kinase (MAPKKK). Typical 
MAPK substrates are cytoplasmic or nuclear proteins, 
such as transcription factors [3]. MAPKKKs are the 
most divergent of these three types of kinases in plants; 
the Arabidopsis thaliana genome contains approxi- 
mately 60 MAPKKKs, 10 MAPKKs, and 20 MAPKs [4]. 

Based on phylogenetic analysis of the amino acid 
sequences of their catalytic kinase domains, plant 
MAPKKKs have been classified into three groups: A, B, and 
C [4]. Group A contains many MAPKKKs involved in PCD 
and stress and defense responses; e.g., A. thaliana 
AtMEKKl is involved in the signaling pathway of basal 
defense induced by PAMPs [5], and Medicago sativa 
MsOMTKl [6] is involved in that of oxidative stress- 
induced cell death. Group A also includes MAPKKKs that 
have important functions in HR induction. Silencing of the 
genes encoding Nicotiana tabacum NPKl (NtNPKl) and 
Nicotiana benthamiana MAPKKKa (NbMAPKKKa) sup- 
presses the N gene-mediated HR induced by the helicase 
domain of tobacco mosaic virus (TMV) replicase and Pto- 
mediated HR induced by Pseudomonas syringae pv. tomato 
(Pst) effector avrPto, respectively [7,8]. Recently, N. 



benthamiana NbMAPKKKe and its tomato {Solanum lyco- 
persicum) ortholog SlMAPKKKe have been implicated in 
PCD induction in the HR against Gram-negative bacterial 
pathogens [9]. In addition, silencing of the genes encoding 
the MAPKK MEK2 and the MAPK SIPK, both of which 
act downstream of NbMAPKKKa, also attenuates the N 
gene-mediated HR against TMV [10]. Conversely, silencing 
of the tomato orthologs of MAPKK MEKl and MAPK 
NTF6, both of whose tobacco orthologs act downstream of 
NtNPKlin tobacco, leads to loss of the Pto-mediated HR in 
tomato [11]. Therefore, the NtNPKl- and NbMAPKKKa- 
initiated MAPK cascades are essential for both the N gene- 
mediated and the Pto-mediated HR, suggesting that at least 
two distinct MAPK cascades are involved in the regulation 
of a single HR event [8]. Furthermore, it is now becoming 
apparent that two distinct MAPK cascades are involved in 
non-HR environmental responses [12]. 

Plants generally appear to use the same MAPKK/ 
MAPK sets in different responses to environmental 
stimuli. The A. thaliana MAPKKs AtMKK4 and AtMKK5 
and/or their downstream component MPK6 are involved 
not only in the signaling pathway for basal defense down- 
stream of AtMEKKl but also in ethylene production and 
stomata formation [5,13,14]. Given the relatively limited 
number of MAPKKs and MAPKs in plants, the diversity 
of these responses (functions) is assumed to be possible 
due to the great diversity of MAPKKKs [15,16]. Therefore, 
comparative functional analysis among MAPKKKs is 
needed to reveal the molecular mechanisms underlying a 
variety of responses to environmental stresses. 

We previously showed that systemic necrosis, the dis- 
ease symptom caused by plantago asiatica mosaic virus 
Lil (PIAMV-Lil), was accompanied by resistance traits 
similar to HR. Using tobacco rattle virus (TRV)-based 
virus-induced gene silencing (VIGS) [17], we demon- 
strated that NbSGTl and NbRARl, which are important 
in the HR, and the MAPK cascade including NbMAPKKKa/ 
NbMEK2, are essential for the induction of PCD-associated 
systemic necrosis induced by PLAMV-Lil [18,19]. This result 
and those described above led us to hypothesize that other 
MAPKKK genes in addition to NbMAPKKKa are involved 
in the systemic necrosis induced by PLAMV-Lil. 

In the present study, we isolated three novel group A 
MAPKKK genes from N. benthamiana, a model plant of 
the family Solanaceae, using an expressed sequence-tag 
(EST) database. The three cloned genes were designated 
NbMAPKKKji, NbMAPKKKy, and NbMAPKKKc2. Fur- 
ther study revealed that NbMAPKKK|3 and NbMAPKKKy 
are positive regulators of PCD. In addition, the results of 
epistasis analysis performed using VIGS and agroinfiltra- 
tion suggest that two of these MAPKKKs (NbMAPKKK|3 
and NbMAPKKKy), together with NbMAPKKKa, com- 
prise a linear signaling pathway important in the induction 
of PCD. 
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Results 

Cloning of three novel group A MAPKKK genes from 
Nicotiana benthamiana 

To conduct a comparative analysis of the roles of 
MAPKKK genes in defense responses, we first sought to 
clone N. benthamiana MAPKKK genes belonging to group 
A. Group A includes four subgroups, A1-A4, and contains 
many genes involved in plant defense responses. To obtain 
one N. benthamiana MAPKKK gene homolog from each 
group A subgroup, we selected the following: the Arabidop- 
sis MAPKKK genes AtMEKKl [20], AtMAPKKKy, and 
AtMAPKKKel [21,22] as representatives of the Al, A2, and 
A4 subgroups, respectively. For the A3 subgroup, which 
includes NtNPKl [23] and AtANPl, we were able to amp- 
lify N. benthamiana cDNA fragments using NtNPKl-spe- 
cific primers; the amplified gene turned out to be 98.0% 
identical in nucleotide sequence to that encoding the 
NtNPKl kinase domain, which has been well characterized 
for its role in defense responses [7]. Therefore, the A3 sub- 
group was excluded from further study. 

Using the nucleotide sequences encoding the highly 
conserved kinase domains of AtMEKKl [19], AtMAPKKKy, 
and AtMAPKKKel as queries for BLAST searches against 
the N. benthamiana, N. tabacum, and Solanum lycopersicum 
EST databases, we obtained three EST sequences: N. 
benthamiana TC15397, N. tabacum BP133312, and S. lyco- 
persicum BI931567, respectively. Based on these sequences, 
we designed specific primers to isolate full-length N. 
benthamiana MAPKKK cDNA clones, as described in the 
Materials and Methods. The cloned N. benthamiana 
MAPKKK genes are predicted to encode proteins with highly 
conserved kinase domains and more divergence in other 
regions, a general feature of plant MAPKKK genes [4]. 

Next, we used the amino acid sequences of the kinase 
domains of our newly cloned MAPKKKs and those of 
other previously identified MAPKKKs to construct a 
phylogenetic tree (Figure lA). Our newly cloned MAPKKK 
genes were tentatively designated after the name of the 
most closely related genes in the phylogenetic tree. Thus, 
our A2 MAPKKK homolog was designated NbMAPKKKy 
after A. thaliana AtMAPKKKy. At a late phase of this 
study, an A4 MAPKKK homolog, N. benthamiana 
NbMAPKKKc, was reported [9]. Although our A4 homolog 
shares high sequence similarity with NbMAPKKKc, the 
genes were not identical (96.6% or 95.0% identity at the nu- 
cleotide or amino acid level, respectively) and differed in 
length. Because these results indicated that A4 MAPKKK 
homologs in N. benthamiana consisted of at least two 
genes, our A4 homolog was designated NbMAPKKKc2 
after NbMAPKKKc. In fact, from the recently released N. 
benthamiana draft genome sequence (http://solgenomics. 
net/), only two contigs that show high sequence homology 
(about -90%) with NbMAPKKKc2 were retrieved and each 
of these sequences corresponded to NbMAPKKKc2 and 



NbMAPKKKc, respectively. Although our Al MAPKKK 
homolog was most closely related to MsOMTKl, which 
was designated after the functional features of its gene pro- 
ducts, designating this gene "NbOMTKl" was apparently 
incorrect, as we do not know if our Al homologous gene 
product and MsOMTKl have a similar function. The next 
most closely related genes were AtMEKKl and Bmssica 
napus BnMAPKKKpi. Therefore, our Al homolog was 
designated NbMAPKKKji after BnMAPKKK^l, a gene of 
the Al subgroup, in accordance with other newly cloned 
genes NbMAPKKKy and NbMAPKKKc2. The NbMAPKKKji, 
NbMAPKKKy, and NbMAPKKKc2 cDNA sequences deter- 
mined in this study were deposited in the DNA Data Bank of 
Japan (DDBJ) under the accession numbers AB649283, 
AB649284, and AB649285, respectively. 

The domain structures of the three novel MAPKKK 
genes are shown in Figure IB. The NbMAPKKKfi, 
NbMAPKKKy and NbMAPKKKc2 cDNA are 2477, 
2794, and 4708 bp in length, respectively, and are re- 
spectively predicted to encode proteins of 564, 651, and 
1395 amino acids. Kinase domains are found in the C- 
terminal regions of NbMAPKKK|3 and NbMAPKKKy 
and the N-terminal region of NbMAPKKKe2. The C- 
terminal region of NbMAPKKKe2 contains two ARM 
(armadillo/|3-catenin-like repeat) domains. The domain 
structures of these proteins are similar to those of their 
respective A. thaliana homologs. 

Overexpression of the NbMAPKKKP and NblVlARKKKy 
kinase domains, but not the NbMAPKKK£2 kinase domain, 
causes cell death 

Overexpression of the kinase domain of the tomato 
homolog of NbMAPKKKa, SlMAPKKKa, in N. benthami- 
ana leaves by agroinfiltration has been reported to induce 
pathogen- independent cell death [8]. To investigate 
whether our three newly identified MAPKKKs also pos- 
sessed cell death-inducing activity, we performed transient 
overexpression analysis of each kinase domain using 
agroinfiltration. The transient overexpression of the kinase 
domains of NbMAPKKKP and NbMAPKKKy induced 
pathogen-independent cell death in the infiltrated area 
(Figure 2A). In the infiltrated areas, cell death was asso- 
ciated with significant hydrogen peroxide production, 
which is detected by a characteristic brown color that 
emerges upon 3,3'-diaminobenzidine (DAB) staining. 
These results were confirmed by measuring the extent of 
cell death using an ion leakage assay (Figure 2B), which 
showed significantly increased ion leakage in the 
NbMAPKKK(3- and NbMAPKKKy-overexpressing areas. 
No significant difference in the level of ion leakage was 
observed for NbMAPKKKP vs. NbMAPKKKy. In contrast, 
overexpression of the kinase domain of NbMAPKKKe2 
failed to induce cell death or hydrogen peroxide produc- 
tion and failed to increase the level of ion leakage. 
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Figure 1 Phylogenetic analysis of plant MAPKKK genes and structural features of novel N. bentbamiana IVIAPKKK genes. A) Phylogenetic 
tree for the kinase domain-encoding regions of group A MAPKKK genes. Numbers represent bootstrap scores. Names of genes used in this study 
are shown in boldface. Names of the three novel IVIAPKKK genes are underlined. B) Diagram of the domain structures of the three novel MAPKKK 
genes. ORFs are indicated by boxes. The 5'- and 3'-untranslated regions (UTRs) are shown as horizontal lines at the left and right, respectively, of 
the boxes. Shaded boxes and gray boxes indicate l<inase domains and ARM domains, respectively. ORFs and UTRs are drawn to the same scale. 



Western blot analysis using anti-Myc monoclonal anti- 
body (Millipore, BUlerica, MA, USA) was performed to as- 
sess the expression levels of these kinase domains 
(Figure 2C). Specific signals of NbMAPKKKp and 
NbMAPKKKe2 kinase domains were detected, but no ac- 
cumulation of NbMAPKKK/ kinase domain was found. 
However, we suppose that the NbMAPKKK/ kinase do- 
main was accurately expressed, even at very low levels, be- 
cause its overexpression rapidly induced cell death. 
Western blot analysis also showed that the accumulated 
expression level of NbMAPKKKe2 kinase domain, whose 
overexpression did not induce cell death, was much 
greater than the expression levels of the NbMAPKKK(3 
and NbMAPKKKy kinase domains. Because these results 
indicated that NbMAPKKKe2 is not involved in the in- 
duction of pathogen-independent cell death, we excluded 
NbMAPKKKe2 from further analysis. 



Overexpression of full-length, catalytlcally active 
NbMAPKKKp or NbMAPKKKy causes cell death 

To further investigate the involvement of NbMAPKKKp 
and NbMAPKKKy in pathogen-independent cell death, 
we conducted transient overexpression experiments 
using their full-length open reading frames (ORFs). As 
shown in Figure 2D, overexpression of NbMAPKKKp or 
NbMAPKKK/ induced cell death. The observed level of 
ion leakage was similar for overexpression of NbMAPKKKp, 
NbMAPKKKy, and NbMAPKKKa, which was used as a 
positive control (Figure 2B). 

Several previous studies have demonstrated that kinase 
activity is required for the triggering of pathogen- 
independent cell death by transient overexpression of a 
component of a MAPK cascade [8,24,25]. To examine 
whether NbMAPKKKp- and NbMAPKKKy-induced ceU 
death also required kinase activity, we constructed fiall- 
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Figure 2 (See legend on next page.) 
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Figure 2 Overexpression of NblVIAPKKK|3 and NbMAPKKKy induces pathogen-independent cell death in N. bentbamiana leaves. N. 

benthamiana leaves were infiltrated witli Agrobacterium strains carrying pEarleyGate203 vector [32] derivatives harboring ihe coding sequences 
for mutant or full-length MAPKKKs. All experiments, except that shown in Figure 4, were repeated at least three times with similar results. A) 
Symptoms of infiltrated W. benthamiana leaf areas overexpressing the indicated MAPKKK kinase domain (KDs) or ATP-binding site-deficient 
(K^ M) KD mutants. Images of the same leaves after DAB staining are shown in the lower panels. DAB staining to detect hydrogen peroxide 
production was performed as previously described [18]. Pictures were taken 5 days post-infiltration (dpi). Each protein was transiently co- 
expressed with silencing suppressor pi 9. Agrobacterium cultures were grown to a turbidity (ODgoo) of 0.5 for use in agroinfiltration. B) Ion leakage 
in the infiltrated leaf areas overexpressing wild-type or M mutant NbMAPKKKf? or NbMAPKKKy or their wild-type or K— > M mutant KDs or 
NbMAPKKK£2 KD. GUS and NbMAPKKKa were used as internal controls. The ion leakage assay was performed as previously described [18]. Data 
shown represent means ± standard deviation of at least three independent plants. C) Western blot analysis of Myc-tagged NbMAPKKKf? or 
NbMAPKKKy KDs or their K^ M mutants or NbMAPKKKe2 KD. Myc-tagged GUS was used as an internal control. Total proteins were extracted 
from each gene-infiltrated area at 5 dpi. Two replicates are shown for each Myc-tagged construct. D) Symptoms of infiltrated leaf areas 
overexpressing full-length wild-type or K^M mutant NbMAPKKKf? or NbMAPKKKy. Images of the same leaves after DAB staining are shown in 
the lower panels. 



length and kinase domain MAPKKK mutants deficient in 
ATP binding. In these mutants, the essential conserved ly- 
sine residue (K) in the ATP-binding site of the kinase do- 
main was replaced with a methionine (M) [8]. As shown in 
Figure 2 A and 2D, none of these K ^ M mutants induced 
cell death when overexpressed. This result was confirmed 
by DAB staining (Figure 2A and 2D) and ion leakage assays 
(Figure 2B). Also, the expression of these kinase domain K 
—^ M mutants was confirmed by western blot analysis with 
anti-Myc monoclonal antibody (Figure 2C). These results 
indicate that cell death induced by the transient overexpres- 
sion of N. benthamiana NbMAPKKK|3 or NbMAPKKKy is 
dependent on their kinase activities. 

Silencing of NbMAPKKKf or NbMAPKKKy suppresses 
PIAMV-Lii -induced PCD 

In light of the above finding that overexpression of full- 
length NbMAPKKKP or NbMAPKKKy or their kinase 
domains can induce cell death, we expected that silen- 
cing of the genes encoding these MAPKKKs would sup- 
press virally induced PCD. Therefore, we used a TRV- 
based VIGS system [17] to silence NbMAPKKKf and 
NbMAPKKKy. Successful silencing of each gene was 
confirmed by analyzing the expression of NbMAPKKKf^, 
NbMAPKKKy, and NbMAPKKKa (control) using real- 
time reverse transcription (RT)-PCR analysis (Figure 3A). 
Also, Southern blot analyses using a kinase domain- 
specific probe of each gene revealed that NbMAPKKKf^ 
and NbMAPKKKy exist as single-copy and multiple-copy 
genes, respectively, in the N. benthamiana genome 
(Additional file 1: Figure SI). In the N. benthamiana draft 
genome, two contigs showing high sequence homology 
(about >90%) with NbMAPKKKy were obtained, suggest- 
ing that there are at least two copies of NbMAPKKKy 
homologs. Thus, we assume that a single gene is specifically 
targeted in NbMAPKKKji-sAenceA plants, but one or more 
highly similar genes might be targeted in NbMAPKKKy- 
silenced plants. Hereafter, we use the term "NbMAPKKKy" 
to indicate NbMAPKKKy and/or its homologous gene(s) 



except for transient overexpression experiments. The gene- 
silenced plants did not display any obvious phenotypic dif- 
ferences compared with non-silenced (VIGS vector) control 
plants (Figure 3B), suggesting that neither NbMAPKKKp 
nor NbMAPKKKy is involved in normal plant growth and 
development. In contrast, silencing of NtNPKl, a 
MAPKKK involved in cell plate formation in plant cytokin- 
esis and in N gene-mediated HR cell death, causes severe 
stunting of plants [7]. 

NbMAPKKKf- and NbMAPKKKy-s\\.encedi plants and 
control non-silenced plants were subsequently inoculated 
with PL\MV-Lil, which induces PCD-associated necrosis 
in N. benthamiana. In the PlAMV-Lil -infiltrated areas of 
non-silenced plants, we observed characteristic necrotic 
symptoms and brown color staining, indicating the accu- 
mulation of hydrogen peroxide. However, in the Lil- 
infiltrated A/fcMAP/<7<7</S-silenced plants, cell death and 
hydrogen peroxide production were completely compro- 
mised. Also, in NbMAPKKKy-sAenceA plants, cell death 
and hydrogen peroxide production were suppressed com- 
pletely and partially, respectively, despite only about 20% 
decrease in the abundance of NbMAPKKKy transcripts 
(Figure 3C). Suppression of cell death in the NbMAPKKK^- 
and NbMAPKKKy-S)\s\\ced plants was confirmed using ion 
leakage assays (Figure 3D). The suppression of cell death in 
A/fcMAP/<7<7<y-silenced plants despite the slight decrease in 
its mRNA can be explained by a strict requirement of a 
high level of NbMAPKKKy protein in cell death. Although 
it is possible that residual NbMAPKKKy protein contri- 
butes to the decreased level of hydrogen peroxide produc- 
tion, its level might be too low to induce cell death. These 
results indicate that NbMAPKKK|3 and NbMAPKKKy are 
involved in hydrogen peroxide production and PCD 
induced by PlAMV-Lil. 

NbMAPKKKa, NbMAPKKKp, and NbMAPKKKp form a 
linear signaling pathway that induces cell death 

Together with the findings of our previous study show- 
ing that NbMAPKKKa is involved in PlAMV-Lil- 
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Figure 3 Silencing of NbMAPKKKP and NbMAPKKKy suppresses PIAMV-Lil-induced PCD. A) Confirmation of specific MAPKKK gene silencing 
in gene-silenced plants. The relative quantities of NbMAPKKKa, NbMAPKKKji, and NbMAPKKKy mRNA in non-silenced or NbMAPKKKa-, NbMAPKKKfi-, 
or Wfa/M/^PKKKy-silenced plants were measured using real-time RT-PCR. The data for each leaf tissue sample were normalized to the Nbl8S rRNA 
gene expression level in the same sample. Data shown represent means ± standard deviation of at least three independent plants. Asterisks 
indicate statistically significant differences from control plants (P<0.05). B) Typical phenotypes observed in non-silenced and NbMAPKKKji- and 
Wfa/M/\P/</(/<y-silenced plants. Pictures were taken 21 dpi. C) Symptoms in the PIAMV-Lil -infiltrated areas of non-silenced and NbMAPKKKfi- or 
NbMAPKKKy-silenced plants (upper panels) and DAB staining of the same leaves (lower panels). Red circles indicate cell death; yellow circles 
indicate no symptoms. Leaves were infiltrated with a PIAMV-Lil -expressing Agrobacterium culture grown to a turbidity (ODgoo) of 0.05. D) Ion 
leakage in the PIAIVIV-Lil -infiltrated areas of non-silenced and NbMAPKKKj}- or NbMAPKKKy- silenced plants. Data shown represent 
means ± standard deviation of at least three independent plants. Asterisks indicate statistically significant differences from control plants (P < 0.05). 



induced PCD [19], the above findings demonstrate that 
NbMAPKKK^ and NbMAPKKKy are also essential for 
this virally induced PCD. To determine the relationships 
among these three genes in the PCD signaling pathway, 
we designed epistasis experiments combining transient 
overexpression and silencing of combinations of these 
genes. In these experiments, plants with a specific VIGS- 
silenced MAPKKK gene were agroinfiltrated with an 
Agrobacterium culture expressing a different MAPKKK 
gene, and the level of ion leakage was measured. Prior to 
the epistasis experiments, a preliminary experiment was 
performed to determine the minimum turbidity of Agro- 
bacterium inoculum for each gene sufficient to induce 
complete cell death. Complete cell death was observed 



at the following cell densities at ODeoo: 0.5, NbMAPKKKa; 
1.0, NbMAPKKKP; and 0.05, NbMAPKKKy. Although each 
MAPKKK was overexpressed at different concentrations of 
Agrobacterium inocula, no significant difference in the basal 
ion leakage level was observed among the different turbid- 
ities (ranging from 0.05 to 1.0) of Agrobacterium inocula 
that expressed the GUS gene in wild-type plants 
(Additional file 2: Figure S2). For NbMAPKKKp, cell death 
induced by both the full-length ORF and kinase domain 
was slow and weak, but cell death induced by the kinase 
domain was more often observed than that induced by the 
full-length ORF. In the ion leakage assay shown in 
Figure 2B, this tendency was detected, although the differ- 
ence in cell death induced by the full-length ORF and 
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kinase domain of NbMAPKKKP was not statistically signifi- 
cant. Previous reports have shown that the substrate specifi- 
city of MAPKKK is not affected when only the kinase 
domain is transiently activated [5,8,26]. Hence, the kinase 
domain of NbMAPKKKp was used for further analysis, in- 
stead of the full-length ORF. 

As shown in Figure 4A, ion leakage levels induced by 
transient expression of the NbMAPKKK|3 kinase domain 
or full-length NbMAPKKKy were lower in NbMAPKKKa- 
silenced plants than in VIGS vector-only control plants. 
Cell death induced by both NbMAPKKKp and 
NbMAPKKKy was consistently significantly suppressed in 
NbMAPKKKa-silenceA plants compared to control plants 
(data not shown). This result suggests that NbMAPKKKa 
either acts downstream of both NbMAPKKKp and 
NbMAPKKKy in cell death induction or plays an essential 
role in codependent activation of both NbMAPKKKp and 
NbMAPKKKy. In contrast, the ion leakage levels induced 
by transient expression of full-length NbMAPKKKa or 
NbMAPKKKy were similar in MMAP/<7<7<)8-silenced 
plants and control plants (Figure 4B), suggesting that 
NbMAPKKKa and NbMAPKKKp are not codependent in 
their activation; rather, NbMAPKKKp functions upstream 
of NbMAPKKKa in the cell death signaling pathway. This 
result led us to postulate that NbMAPKKKp functions up- 
stream of NbMAPKKKy. As expected, in NbMAPKKKy- 
silenced plants, cell death and the ion leakage level 
induced by the transient overexpression of the 
NbMAPKKKp kinase domain were significantly compro- 
mised compared to control plants (data not shown and 
Figure 4C, respectively). Cell death and the ion leakage 
induced by transient overexpression of NbMAPKKKa, 
however, were similar in iVfeAL4P/<7<7<y-silenced and con- 
trol plants (data not shown and Figure 4C, respectively), 
confirming that NbMAPKKKp functions upstream of 
NbMAPKKKy and that NbMAPKKKy functions upstream 
of NbMAPKKKa. In this epistasis analysis, ion leakage 
levels were closely associated with the intensity of cell 
death at all combinations of these genes. Thus, these three 
MAPKKK genes form a linear signaling pathway leading 
to PCD in which NbMAPKKKp and NbMAPKKKa func- 
tion as the furthest upstream and downstream compo- 
nents, respectively, of the three MAPKKK components. 

Discussion 

In our previous study [19], we showed that a MAPK cas- 
cade including NbMAPKKKa and NbMEK2 is required 
for both the PCD-associated systemic necrosis induced 
by PIAMV-Lil and the 7?;«:-mediated HR against potato 
virus X. Several studies have demonstrated that at least 
two MAPK cascades are involved in plant responses to 
various environmental stimuli [8,12]. The great diversity 
of plant MAPKKK genes is assumed to underlie the abil- 
ity of plants to mount specific signaling responses to 



various environmental stimuli [15,16]. Therefore, to bet- 
ter understand the functions of MAPKKKs in the induc- 
tion of PCD, we performed a comparative functional 
analysis of N. benthamiana MAPKKK genes. In the pre- 
vious and present study, we demonstrated that three of 
these genes, NbMAPKKKa [19], NbMAPKKKp, and 
NbMAPKKKy (and/or its homologs), function as positive 
regulators of PlAMV-Lil-induced PCD and are not 
functionally redundant in the cell death signaling path- 
way. In the VIGS experiments, not only NbMAPKKKy 
but also its homologous gene(s) were presumably 
knocked down in TWV-NbMAPKKKy-iniecteA plants. 
However, only one gene in the A2 subgroup was 
obtained from the cloning experiments, suggesting that 
NbMAPKKKy homologous gene(s), which carries almost 
the same sequence as NbMAPKKKy itself, must be ef- 
fectively knocked down in TRV-7V^&MAP/<7<7<r/-infected 
plants. In addition, overexpression of the NbMAPKKKy 
full-length ORF and kinase domain could induce rapid 
cell death. Therefore, NbMAPKKKy itself could play an 
important role in the cell death signaling pathway. 

We also explored the functional relationship among 
these three MAPKKK genes by performing an epistasis 
experiment based on the assumption that if signaling 
protein A functions upstream of signaling protein B, sig- 
naling by activated A will be suppressed by silencing of 
B expression, whereas signaling by activated B will not 
be suppressed by silencing of A expression. Our results 
(Figure 4) demonstrate that NbMAPKKKa, NbMAPKKKp, 
and NbMAPKKKy, which represents NbMAPKKKy and/or 
its homologous gene(s), form a linear signaling pathway 
leading to cell death induction in which NbMAPKKKp acts 
upstream of NbMAPKKKy and NbMAPKKKy acts up- 
stream of NbMAPKKKa. Results of a similar epistasis ana- 
lysis have suggested that an NPK1/MEK1/NTF6 MAPK 
cascade functions downstream of a MAPKKKa/MEK2/ 
SIPK MAPK cascade in the induction of PCD in plants [8]. 
These results support the possibility that NbMAPKKKa, 
NbMAPKKKp, NbMAPKKKy, and NPKl initiate four dis- 
tinct MAPK cascades that are coordinately involved in 
plant cell death. Alternatively, different sets of MAPK cas- 
cades might be involved in cell death signaling pathways 
induced by different plant-microbe interactions. To investi- 
gate this issue, further research is needed to determine 
whether NbMAPKKKp and NbMAPKKKy participate in 
the N gene-mediated HR against TMV or in the Pto- 
mediated HR against Pst. 

In our experiments, the transient overexpression of 
the NbMAPKKKe2 kinase domain by agroinfiltration 
did not induce cell death (Figure 2A). However, Melech- 
Bonfil and associates have shown that the tomato homo- 
log of NbMAPKKKe2, SlMAPKKKe, is required for PCD 
induction in the HR against Gram-negative bacterial 
pathogens [9]. Furthermore, they showed that the 
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Figure 4 Epistasis analysis of the functional relationships 
among NbMAPKKKa, NblVlAPKKKp, and NblMAPKKKv. Epistasis 
analysis was performed by combining functional activation and 
suppression of specific combinations of MAPKKKs using transient 
overexpression and VIGS, respectively. Cell death in infiltrated areas 
was quantified using ion leakage assays at 4 dpi. Data shown 
represent means ± standard deviation for at least five plants. 
Asterisks indicate statistically significant differences from control 
plants (P<0.05). All experiments described in this figure were 
repeated two times with similar results. A) Cell death in the 
NbMAPKKKp kinase domain (pkinase)- or NbMAPKKKy (yORF)- 
overexpressing areas of non-silenced and NbMAPKKKa-sWenced 
plants. Agrobactehum cultures expressing pkinase and yORF were 
grown to a turbidity of 1.0 and 0.05, respectively, for infiltration. B) 
Cell death in the NbMAPKKKa (aORF)- or yORF-overexpressing areas 
of non-silenced and NhMARKKKfi-sWericed plants. Agrobacterium 
cultures expressing aORF and yORF were grown to a turbidity of 0.5 
or 0.05, respectively, for infiltration. C) Cell death in the aORF- or 
pkinase-overexpressing areas of non-silenced and NbMAPKKKy- 
silenced plants. Agrobacterium cultures expressing aORF and pkinase 
were grown to a turbidity of 0.5 or 1.0, respectively, for infiltration. 



transient overexpression of the tomato SlMAPKKKe kin- 
ase domain in N. benthamiana induces pathogen- 
independent cell death. A comparison between the amino 
acid sequences of the kinase domains of NbMAPKKKe2 
and SlMAPKKKe indicated only three amino acid differ- 
ences. These three residues are not highly conserved and 
are not predicted to be essential for plant kinase catalytic 
activity, so explaining this functional difference solely by the 
amino acid sequence level is difficult. Although at first 
glance, our results seem to contradict theirs, the results 
cannot be direcdy compared because three factors in our 
transient overexpression experiments differed from theirs: 
promoters, Agrobacterium strains, and species origin of the 
homologous genes. These three factors, which could influ- 
ence the expression levels of overexpressed proteins, are 
critical for the execution of cell death elicited by proteins 
possessing cell death-inducing activity [27,28]. Moreover, 
the amount of overexpressed NbMAPKKKe2 kinase do- 
main in our experiment was much more than the amounts 
of NbMAPKKKp kinase domain and its K ^ M mutant or 
NbMAPKKKy kinase domain that can induce cell death 
and was comparable to the amount of the NbMAPKKKy 
kinase domain K ^ M mutant. Therefore, NbMAPKKKe2 
might be a less potent inducer of cell death than other 
MAPKKKs. 

Our suggestion that several MAPK cascades function 
in a linear signaling pathway in the induction of cell 
death raises a concern about their roles in plant defense 
responses. To prevent pathogen invasion, plants com- 
monly employ two types of receptors [NBS-LRR (nu- 
cleotide-binding site and leucine-rich repeat) receptors 
and RLK (receptor-like kinase)-type receptors] [1] to 
sense multiple PAMPs or effectors derived from bio- 
trophic and necrotrophic pathogens. When plants 
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perceive pathogens, the transcription of a common set 
of genes that act against various types of pathogens is 
activated [29,30], subsequently producing many defense 
responses, including cell death and the production of 
low-molecular-weight signaling compounds such as 
ethylene and reactive oxygen species (ROS). Activation 
of the MAPKKKa/MEK2/SIPK cascade can induce pro- 
duction of ethylene and ROS [31,32]. NPKl, which initi- 
ates another MAPK cascade, is activated by the ROS 
produced by the MAPK cascade initiated by MAPKKKa 
[26]. Similarly, the alfalfa MAPKKK MsOMTKl is acti- 
vated by hydrogen peroxide [6]. These results suggest 
that compounds induced by an upstream-acting MAPK 
cascade function as signaling molecules that activate a 
downstream-acting MAPK cascade. In addition, they 
suggest that the involvement of several MAPK cascades 
in defense responses enables plants to activate various 
responses simultaneously or coordinately to combat 
various types of pathogens. 

Conclusions 

In this study, we demonstrated that three MAPKKKs in 
N. benthamiana form a linear signaling pathway leading 
to PCD, implying that the involvement of multiple 
MAPK cascades in plant defense responses enables 
plants to exhibit various reactions simultaneously. Fur- 
ther analysis to identify downstream MAPKK and 
MAPK genes directly phosphorylated by NbMAPKKK(3 
and NbMAPKKKy should help clarify the regulatory 
mechanisms of cell death involving these MAPKKK 
genes. Furthermore, to elucidate the mechanisms under- 
lying the sequential activation of each MAPK cascade 
during plant defense responses, the upstream compo- 
nents or chemical compounds that directly activate these 
MAPKKKs must be identified. 

Methods 

Plant materials and virus isolate 

N. benthamiana plants were grown in a growth chamber 
at 25 °C. To virally induce PCD, plants were inoculated 
with the binary vector pLil, which contains the full- 
length cDNA of the Lil isolate of plantago asiatica mo- 
saic virus (PIAMV) downstream of the cauliflower mo- 
saic virus 35 S promoter [33]. 

Cloning of MAPKKK genes from N. benthamiana 

To retrieve partial MAPKKK cDNA sequences, a BLAST 
search was performed using the amino acid sequences of 
the kinase domains of the Arabidopsis MAPKKK genes 
AtMEKKl, AtMAPKKKy, and AtMAPKKKel as queries 
against the N. benthamiana, N. tabacum, and S. lycopersi- 
cum databases of the Gene Index Project at the Computa- 
tional Biology and Functional Genomics Laboratory Web 
site (http:/ / compbio.dfd.harvard.edu/ tgi/ cgi-bin/ tgi/Blast/ 



index.cgi) and the tobacco BY-2 EST clone database of the 
RIKEN BioResource Center (http://www.brc.riken.go.jp/ 
lab/epd/blast/index.shtml). The consensus nucleotide 
sequences for the retrieved partial cDNA sequences and 
query sequences were used to design a pair of specific pri- 
mers for each of the three MAPKKK genes. Each pair of 
specific primers was used for RT-PCR amplification of a 
cDNA fragment from total RNA extracted from N. 
benthamiana leaves. The RT-PCR products were gel- 
purified and cloned into the pGEM-T easy vector (Pro- 
mega, Madison, WI, USA) for sequencing. New primers 
were then designed based on the regions of the obtained 
cDNA sequences that were nonhomologous to the query 
sequences (to prevent nonspecific amplification) and used 
for specific 5 - and 3 -rapid amplification of cDNA ends 
(RACE) using a GeneRacer Kit (Invitrogen, Carlsbad, CA, 
USA) according to the standard protocol provided by the 
manufacturer. The sequences of primers used in this study 
are listed in Table L cDNA sequences of these MAPKKK 
genes were determined from at least three independent 
clones. 

ORE and motif analyses were performed using the ORF 
Finder (http://www.ncbi.nlm.nih.gov/gorf/gorfhtml) and 
CD-Search (http://www.ncbi.nlm.nih.gov/ Structure/ cdd/ 
wrpsb.cgi) programs, respectively. Phylogenetic analysis 
was performed using MEGA3.1 based on a multiple align- 
ment created using ClustalW. The MAPKKK gene 
sequences used in phylogenetic analysis are listed in the 
GenBank database under the following accession num- 
bers: AtMEKKl (NM_1 16919), AtMEKK2 (NM_1 16917), 
AtMEKKS (NM_1 16916), AtMEKK4 (NM_1 17272), 
BnMAPKKKfil (AJ010093), AtANPl (NM_100771), 
AtANP2 (NM_104370), AtANPS (NM_1 11477), NtNPKl 
(D26601), AtMAPKKKy (NM_126084), AtYODA 
(AY357949), NbMAPKKKa (AY500155), AtMAPKKKa 
(NM_179472), AtMAPKKKel (NM_1 12199), AtMAPKKKc2 
(NM_111677), MsOMTKl (AJ575100), BnMAPKKKal 
(AJ010091), BnMAPKKKcl (AJ238845), SlMAPKKKc 
(GU192457), and NbMAPKKKc (GU205153). Homo sapiens 
HsMEKKl (AF042838) was used as an outgroup. 

Construction of plasmids for transient expression of wild- 
type and mutant MAPKKKs and MAPKKK kinase domains 
and agroinfiltration 

The coding regions of NbMAPKKKa, NbMAPKKK^, 
and NbMAPKKKy were amplified by RT-PCR using 
primer pairs Kp-NbMAPKKKa-lF/Xh-NbMAPKKKa- 
1836R, Kp-NbMAPKKKb-lF/EcV-NbMAPKKKb-1695R, 
and Kp-NbMAPKKKg-lF/Xh-NbMAPKKKg-1956R, re- 
spectively. The kinase domains of NbMAPKKK/i, Nb 
MAPKKKy, and NbMAPKKKc2 were amplified using 
primer pairs Kp-NbMAPKKKbkinase-F/EcV-NbMAPKK 
Kbkinase-R, Kp-NbMAPKKKgkinase-F/Xh-NbMAPKKK 
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Table 1 Primers used in this study 



Name 



Sequence (5'-3') 



Comments 



NbTC9992-1F 

NbTC9992-1280R 

NtBPI 333-1 F 

NtBP1333-447R 

LeBI9315-123F 

LeBI9315-658R 

NbTC9992-458R-5RACE 

NbTC9992-49R-5RACE 

NbTC9992-Fl 

NtMAP3Kb-566F 

NtBPI 333-1 99R-5RACE 

NtBPI 333-91 R-nested 

NtBPI 333-333F-5RACE 

NtBPI 333-375F-nested 

NtBPI 333-F1 

NtBPI 333-R1 

LeBI9315-507R-5RACE 

LeBI9315-155R-nested 

LeBI9315-482F-3RACE 

LeBI9315-564F-nested 

LeBI9315-F1 

LeBI9315-F2 

LeBI9315-F3 

LeBI9315-Rl 

LeBI9315-R2 

Kp-NbMAPKKKa-IF 

Xh-NbMAPKKKa-1836R 

Kp-NbMAPKKKb-IF 

EcV-NbMAPKKKb-1695R 

Kp-NbMAPKKKg-IF 

Xh-NbMAPKKKg-1956R 

Kp-NbMAPKKKbkinase-F 

EcV-NbMAPKKKbkinase-R 

Kp-NbMAPKKKgkinase-F 

Xh-NbMAPKKKgkinase-R 

Kp-NbMAPKKKekinase-F 

Xh-NbMAPKKKekinase-R 

NbMAPKKKb-K381M-F 

NbMAPKKKb-K381M-R 

NbMAPKKKg-K374M-F 

NbMAPKKKg-K374M-R 

NbMAPKKKa-255F 

NbMAPKKKa-393R 

NbMAPKKKb-1179F 



GCTGTCAAAGAAGTGTCATTA 

ACCGrrTArTAATCACTATATTGC 

CTTAATGGGCAAGCAGCTAATC 

TCAAGATTGTATGTTGTCTGCTC 

GTTGCAATTAAACAAGTTTGTCTGGA 

GGGTGAAGATCATAGTAGGG 

GCrrGTCCATAGCCTTGGTTCTTCCT 

TTTGCCTTCCCGCATCGCCTTGAT 

GGArTGGAAAGGGGGAAGCT 

TCCGCCGGTCATGTCACT 

ACATGGCTGCAGCTGGTTCATATTC 

TACTATCCTTCTGCATAACTGACTGCAA 

GAGAATCTGTGATCAGCCAGATGTTC 

GAAACCTGTTGGTGGGGTACGAAT 

AGACGCGCATAATTCGCATC 

TGAGCTCTCGrrTGGTAATAAG 

CATCTGCCTCTGTGAACnTGTTGGA 

GaCCrGAGCAATATTCTCCAGAGA 

TGCAACAAAGrrGACAGAGGCAGATG 

ATGTCGGGAGTATGTGCTGCATCTG 

CTGAGAAAGTTC7TGGAAACG 

GAACAGATGAATCAGAAGATG 

AAAGCACTGCATATAAACACAG 

GGTATACATCAAGTCCACGAT 

GAGGGAGTATGCTCTCATG 

GGGGTACCGAATGCCTGCTTGGTGGGGAA 

GGGTGGAGTGCTAAAGAATTGGTCTTAGTTrTG 

GCGGTAGCGAATGCATCGATTGCCAGGAArmTGC 

GGGATATCTTTAAAGCCTCTTGCCCAGAnTTG 

GGGGTACCGAATGCGTTGGTGGCAGAACG 

GGGTGGAGTGCTACCTCrCTAGAGATAAACG 

GCGGTACCGAATGTAGTGGGACAAAGGTGATCT 

GGGATATCTTTACACAAAAGGATGCTCCAAGA 

GGGGTACCGAATGTGGCAAAAAGGGAAGCTTATTG 

GGCTCGAGTGrTACATAAATGGATGTTCCAATAAG 

GGGGTACCGAATGAAATATATGCrCGGAGATGAG 

GGGTGGAGTGTTATATCGATGGATGTGAAAGTAG 

"mTGCTGTCATGGAAGTGTCATTACTTGATCA 

ATGACACTTCCATGACAGCAAAAAAGAAACCG 

CTGGAGCTTTATGTGCGATGATGGAAGTrGAATT A1TACCGGA 

TCGGGTAATAATTGAACTTCCATCATCGCACATAA AGCTCGAG 

GGTTGirrTGGGATGTGGGGTCAG 

GAGTGGGCTCAACCTATTATCGCC 

GACAAGGCAGATTTTACATGGTTrG 



Specific primer for TC9992 

Specific primer for TC9992 

Specific primer for BP133312 

Specific primer for BP133312 

Specific primer for BI931567 

Specific primer for BI931567 

5'-RACE for TC9992 

5'-RACE for TC9992 

Sequencing 

Sequencing 

5'-RACEfor BP133312 

5'-RACEfor BP133312 

3'-RACEfor BP133312 

3'-RACEfor BP133312 

Sequencing 

Sequencing 

5'-RACEfor LeBI931567 

5'-RACE for LeBI931567 

3'-RACE for LeBI931567 

3'-RACEfor LeBI931567 

Sequencing 

Sequencing 

Sequencing 

Sequencing 

Sequencing 

Full-length ORE of NbMAPKKKa 
Full-length ORE of NbMAPKKKa 
Full-length ORE of NbMAPKKKp 
Full-length ORE of NbMAPKKK(3 
Full-length ORE of MbMARKKKy 
Full-length ORE of MbMARKKKy 
Kinase domain of NbMAPKKKp 
Kinase domain of NbMAPKKKp 
Kinase domain of NbMAPKKKy 
Kinase domain of NbMAPKKKy 
Kinase domain of NbMAPKKK£2 
Kinase domain of NbMAPKKK£2 
K^ M mutant of NbMAPKKKfJ 
K^ M mutant of NbMAPKKKp 
K^ M mutant of NbMAPKKKy 
K^ M mutant of NbMAPKKKy 
Real-time RT-PCR for NbMAPKKKa 
Real-time RT-PCR for NbMAPKKKa 
Real-time RT-PCR for NbMAPKKKp 
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Table 1 Primers used in this study (Continued) 



NbMAPKKKb-1286R 


AGCTTGACCGATCCGTTAGCA 


Real-time RT-PCR for NbMAPKKKp 


NbMAPKKKg-903F 


CCGTGAGTGTAGTGCTCAGGGTAA 


Real-time RT-PCR for NbMAPKKKy 


NbMAPKKKg-984R 


TGCCGTAGGCTGGTGTGATG 


Real-time RT-PCR for NbMAPKKKy 


Nb18S-193F 


ATACGTGCAACAAACCCCGAC 


Real-time RT-PCR for NblBS rRNA 


Nb18S-280R 


TGAATCATCGCAGCAACGG 


Real-time RT-PCR for NblSS rRNA 



gkinase-R, and Kp-NbMAPKKKekinase-F/Xh-NbMAPK 
KKekinase-R, respectively. 

Loss-of-function mutants of NbMAPKKKp and 
NbMAPKKKy were obtained by substituting methionine 
for the essential lysine in the ATP-binding site in the 
kinase domain. These mutants were obtained by site- 
directed mutagenesis using primers NbMAPKKKb- 
K381M-F and NbMAPKKKb-K381M-R for NbMAPKKKp 
and NbMAPKKKg-K374M-F and NbMAPKKKg-K374M- 
R for NbMAPKKKy. Each PCR-amplified MAPKKK gene 
fragment was subcloned into the pEarleyGate 203 vector 
under 35 S promoter [34] via the pENTA entry vector 
[35] with LR Clonase II Enzyme Mix (Invitrogen). GUS 
gene was also subcloned into the pEarleyGate 203 vector 
under the 35 S promoter. These plasmid vectors were 
transformed into Agrobacterium tumefaciens strain 
EHA105. Agroinfiltration was performed as previously 
described [27]. 

Construction of VIGS vectors and VIGS 

VIGS was performed as described previously using 
pTV:00 derivatives and pBintra6 [17]. NbMAPKKKa si- 
lencing was induced using pTViNbMAPKKKa [19]. For 
NbMAPKKKp and NbMAPKKKy silencing, the 390-bp 
KpnI-PvuII fragment of the 5 -terminal region of the 
PCR-amplified full-length NbMAPKKKP ORE and the 
576-bp KpnI-EcoRV fragment of the PCR-amplified 
full-length NbMAPKKKy ORE were introduced into 
pTV:00 in the antisense direction to generate pTV: 
NbMAPKKK^ and pTY-.NhMAPKKKy, respectively. The 
NbMAPKKKp ORE was amplified using primers Kp- 
NbMAPKKKb-1 F and EcV-NbMAPKKKb-1695R, and 
the NbMAPKKKy ORE was amplified using primers Kp- 
NbMAPKKKg-1 E and Xh-NbMAPKKKg-1956R. 

Real-time RT-PCR 

The methods used for total RNA isolation, cDNA synthesis, 
and quantitative analysis of gene expression were the same 
as those previously described [18]. Expression of 18 S rRNA 
was used as the normalizer. The primers used for real- 
time RT-PCR were NbMAPKKKa-255E and NbMAPKK 
Ka-393R for NbMAPKKKa, NbMAPKKKb-1179E and 
NbMAPKKKb-1286R for NbMAPKKK|3, NbMAPKKKg- 
903E and NbMAPKKKg-984R for NbMAPKKKy, and 
Nbl8S-193E and Nb280R for Nbl8S rRNA. 



Additional files 



Additional file 1: Figure SI. Souttiern blot analysis of A) NbMAPKKKfi 
and B) NbMAPKKKy using kinase domain-specific DNA probes. DNA 
probes were generated by using the PCR DIG Probe Synthesis Kit (Roche, 
Basel, Switzerland) according to the manufacturer's instructions. Each lane 
was loaded with 5 |ig of total genomic DNA digested with each 
restriction enzyme. 

Additional file 2: Figure S2. Ion leakage of GUS-infiltrated areas. An 
Agrobacterium strain expressing the GUS gene using the 35 S promoter 
was infiltrated at the following turbidities: 0.05, 0.5, and 1.0. 
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MR: Hypersensitive response; (VIAPK: (Vlitogen-activated protein kinase; 
IVIAPKK: IVIAPK kinase; IVIAPKKK MAPKK kinase; PCD: Programmed cell death; 
PIAMV: plantago asiatica mosaic virus; Pst: Pseudomonas syringae pv. tomato; 
TMV: tobacco mosaic virus. 

Authors' contributions 

IVIH, KK, YY and SN designed the experiments. IVIH and KK performed the 
experiments and analyzed the data together with KIVl, YO, YY and SN. TS, Kl 
and YT contributed new reagents and analytic tools. MH, KK and SN wrote 
the paper All authors discussed the results and approved the final 
manuscript. 

Acl<nowledgements 

We thank Dr David Baulcombe (University of Cambridge, Cambridge, UK) for 
tobacco rattle virus vector and silencing suppressor pi 9 gene. This work 
supported by grand-in-aid from the Japan Society for the Promotion of 
Science and, in part, by a Grand-in-aid for JSPS Fellows, and by the Program 
for Promotion of Basic Research Activities for Innovative Bioscience 
(PROBRAIN). 

Received: 19 IVIarch 2012 Accepted: 26 June 2012 
Published: 8 July 2012 

References 

1. Jones JD, DangI JL: The plant immune system. Nature 2006, 444 
(7117):323-329. 

2. Rogers HJ: Cell death and organ development in plants. Curr Top Dev Biol 
2005, 71:225-261. 

3. Rodriguez IVIC, Petersen IVI, IVlundy J: Mitogen-activated protein kinase 
signaling in plants. Anr^u Rev Plant Biol, 61:621-649. 

4. MAPKGroup: Mitogen-activated protein kinase cascades in plants: a new 
nomenclature. Trends Plant Scl 2002, 7(7):301-308. 

5. Asai T, Tena G, Plotnikova J, Willmann IVIR, Chiu WL, Gomez-Gomez L, Boiler 
T, Ausubel FM, Sheen J: MAP kinase signalling cascade in Arabidopsis 
innate immunity. Nature 2002, 415(6875):977-983. 

6. Nakagami H, Kiegerl S, Hirt H: OIVITKI, a novel MAPKKK, channels 
oxidative stress signaling through direct MARK interaction. J Biol Ciiem 
2004 279(26):26959-26966. 

7. Jin H, Axtell IVU, Dahlbeck D, Ekwenna 0, Zhang S, Staskawicz B, Baker B: 
NRK1, an MEKKI-like mitogen-activated protein kinase kinase kinase, 
regulates innate immunity and development in plants. Dev Cell 2002, 3 
(2):291-297. 

8. del Pozo 0, Pedley KF, Martin GB: MARKKKa is a positive regulator of cell 
death associated with both plant immunity and disease. Et^BO J 2004, 23 
(15)3072-3082. 



Hashimoto ef al. BMC Plant Biology 201 2, 1 2:1 03 Page 1 3 of 1 3 

http://www.biomedcentral.eom/1 471 -2229/1 2/1 03 



10. 



12. 



13. 



14. 



16. 



17. 



19. 



20. 



22. 



23. 



24. 



25. 



26. 



27. 



28. 



29. 



Melech-Bonfil S, Sessa G: Tomato IVIAPKKKe is a positive regulator of cell- 
death signaling networks associated with plant immunity. Plant J 2010, 

64(3):379-391. 

Jin H, Liu Y, Yang KY, Kim CY, Baker B, Zhang S: Function of a mitogen- 
activated protein kinase pathway in N gene-mediated resistance in 
tobacco. Plant J 2003, 33(4):719-731. 

Ekengren SK, Liu Y, Schiff M, Dinesh-Kumar SP, Martin GB: Two MARK 
cascades, NPRl, and TGA transcription factors play a role in Pto- 
mediated disease resistance in tomato. Plant J 2003, 36(6):905-917. 
Yoo SD, Cho YH, Tena G, Xiong Y, Sheen J: Dual control of nuclear EIN3 by 
bifurcate MAPK cascades in C2H4 signalling. Nature 2008, 451 {7180):789-795. 
Liu Y, Zhang S: Phosphorylation of l-aminocyclopropane-1-carboxylic acid 
synthase by MPK6, a stress-responsive mitogen-activated protein kinase, 
induces ethylene biosynthesis in Arabidopsis. Plant Cell 2004, 16(12):3386-3399. 
Wang H, Ngwenyama N, Liu Y, Walker JC, Zhang S: Stomatal development 
and patterning are regulated by environmentally responsive mitogen- 
activated protein kinases in Arabidopsis. Plant Cell 2007, 19(l):63-73. 
Zhang S, Klessig DF: IVIAPK cascades in plant defense signaling. Trends 
Plant Sci 200], 6(1 1):520-527. 

Pedley KF, Martin GB: Role of mitogen-activated protein kinases in plant 
immunity. Curr Opin Plant Biol 2005, 8{5):541-547. 
Ratcliff F, Martin-Hernandez AM, Baulcombe DC: Tobacco rattle virus as a 
vector for analysis of gene function by silencing. Plant J 2001 , 25(2):237-245. 
Komatsu K, Yamaji Y, Ozek J, Hashimoto M, Kagiwada 5, Takahashi S, 
Namba S: Nucleotide sequence analysis of seven Japanese isolates of 
Plantago asiatica mosaic virus (PIAMV): a unique potexvirus with 
significantly high genomic and biological variability within the species. 
Arch Virol 2008, 153(1):193-198. 

Komatsu K, Hashimoto M, Ozeki J, Yamaji Y, Maejima K, Senshu H, Himeno 
M, Okano Y, Kagiwada S, Namba S: Viral-induced systemic necrosis in 
plants involves both programmed cell death and the inhibition of viral 
multiplication, which are regulated by independent pathways. Mol Plant 

Microbe Interact 201 0 23(3):283-293. 

Mizoguchi T, Irie K, Hirayama T, Hayashida N, Yamaguchi-Shinozaki K, 
Matsumoto K, Shinozak K: A gene encoding a mitogen-activated protein 
kinase kinase kinase is induced simultaneously with genes for a 
mitogen-activated protein kinase and an S6 ribosomal protein kinase by 
touch, cold, and water stress in Arabidopsis thaliana. Proc Natl Acad Sci 
USA 1996, 93(2):765-769 

Jouannic S, Champion A, Segui-Simarro JM, Salimova E, Picaud A, Tregear J, 
Testillano P, Risueno MC, Simanis V, Kreis M, ef al: The protein kinases 
AtMAP3Kel and BnMAP3Ke1 are functional homologues of S. pombe 
cdc7p and may be involved in cell division. Plant J 2001 , 26(6):637-649. 
Chaiwongsar S, Otegui MS, Jester PJ, Monson SS, Krysan PJ: The protein 
kinase genes MAP3Ke1 and MAP3Ke2 are required for pollen viability in 
Arabidopsis thaliana. Plant J 2006, 48(2): 193-205. 

Banno H, Hirano K, Nakamura T, irie K, Nomoto S, Matsumoto K, Machida Y: 
NPK1, a tobacco gene that encodes a protein with a domain 
homologous to yeast BCK1, STE11, and Byr2 protein kinases. Mol Cell Biol 
1993, 13(8):4745-4752. 

Ren D, Yang KY Li GJ, Liu Y, Zhang S: Activation of Ntf4, a tobacco 
mitogen-activated protein kinase, during plant defense response and its 
involvement in hypersensitive response-like cell death. Plant Physiol 2006, 
141(4):1482-1493. 

Zhang S, Liu Y: Activation of salicylic acid-induced protein kinase, a 
mitogen-activated protein kinase, induces multiple defense responses in 
tobacco. Plant Cell 2001, 13(8):1 877-1 889. 

Kovtun Y, Chiu WL, Tena G, Sheen J: Functional analysis of oxidative 
stress-activated mitogen-activated protein kinase cascade in plants. Proc 
Natl Acad Sci USA 2000, 97(6):2940-2945. 

Komatsu K, Hashimoto M, Maejima K, Shiraishi T, Neriya Y, Miura C, Minato 
N, Okano Y, Sugawara K, Yamaji Y, Namba S: A necrosis-inducing elicitor 
domain encoded by both symptomatic and asymptomatic Plantago 
asiatica mosaic virus isolates, whose expression is modulated by virus 
replication. Mol Plant Microbe Interact 201 1, 24(4):408-420. 
Zhang Y, Dorey S, Swidersk M, Jones JD: Expression of RPS4 in tobacco 
induces an AvrRps4-independent HR that requires EDS1, SGT1 and 
HSP90. Plant J 2004 40(2):21 3-224. 

Navarro L, Zipfel C, Rowland 0, Keller I, Robatzek S, Boiler J, Jones JD: The 
transcriptional innate immune response to flg22. Interplay and overlap 



30. 



32. 



33. 



34. 



35. 



with Avr gene-dependent defense responses and bacterial 
pathogenesis. Plant Physiol 2004 1 35(2):1 113-11 28. 
Tao Y, Xie Z, Chen W, Glazebrook J, Chang HS, Han B, Zhu T, Zou G, Katagiri 
F: Quantitative nature of Arabidopsis responses during compatible and 
incompatible interactions with the bacterial pathogen Pseudomonas 
syringae. Plant Cell 2003, 15(2):31 7-330 

Kim CY, Liu Y, Thorne ET, Yang H, Eukushige H, Gassmann W, Hildebrand D, 
Sharp RE, Zhang S: Activation of a stress-responsive mitogen-activated 
protein kinase cascade induces the biosynthesis of ethylene in plants. 

Plant Cell 2003, 15(1 1):2707-2718 

Liu Y, Ren D, Pike S, Pallardy S, Gassmann W, Zhang S: Chloroplast- 
generated reactive oxygen species are involved in hypersensitive 
response-like cell death mediated by a mitogen-activated protein kinase 
cascade. Plant J 2007, 51(6):941-954 

Ozek J, Takahashi S, Komatsu K, Kagiwada S, Yamashita K, Mori T, Hirata H, 
Yamaji Y, Ugaki M, Namba S: A single amino acid in the RNA-dependent 
RNA polymerase of Plantago asiatica mosaic virus contributes to 
systemic necrosis. Arch Virol 2006, 151(10):2067-2075. 
Farley KW, Haag JR Pontes 0, Opper K, Juehne T Song K, Pikaard CS: 
Gateway-compatible vectors for plant functional genomics and 
proteomics. Plant J 2006, 45(4):61 6-629. 

Himeno M, Maejima K, Komatsu K, Ozeki J, Hashimoto M, Kagiwada S, 
Yamaji Y, Namba S: Significantly low level of small RNA accumulation 
derived from an encapsidated mycovirus with dsRNA genome. Virology, 
396(1 ):69-75. 



doi:1 0.1 1 86/1 471-2229-1 2-1 03 

Cite this article as: Hashimoto ef a/.: Identification of three MAPKKKs 
forming a linear signaling pathway leading to programmed cell death 
in Nicotiana benthamiana. BMC Plant Biology 20]2 12:103. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



